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ABSTRACT 



Active galaxies are thought to be both fueled and obscured by neutral gas 
removed from the host galaxy and funneled into a central accretion disk. We 
performed a VLA imaging survey of 21 cm absorption in Seyfert and starburst 
nuclei to study the neutral gas in the near- nuclear environment. With the 
exception of NGC 4151, the absorbing gas traces 100 pc scale, rotating disks 
aligned with the outer galaxy disk. These disks appear to be rich in atomic gas 
relative to nuclear disks in non-active spirals. We find no strong evidence for 
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rapid in-fall or out-flow of neutral hydrogen, but our limits on the mass infall 
rates are compatible with that required to feed a Seyfert nucleus. 

Among the galaxies surveyed here, neutral hydrogen absorption traces 
parsec-scale gas only in NGC 4151. Based on the kinematics of the absorption 
line, the disk symmetry axis appears to align with the radio jet axis rather 
than the outer galaxy axis (cf. [Mundell et al. 1995|) . The most surprising 
result is that we detect no 21 cm absorption towards the central radio sources 
of the hidden Seyfert 1 nuclei Mrk 3, Mrk 348, and NGC 1068. Moreover, 
21 cm absorption is commonly observed towards extended radio jet structure 
but appears to avoid central, compact radio sources in Seyfert nuclei. To 
explain these results, we propose that 21 cm absorption towards the nucleus is 
suppressed by either free- free absorption, excitation effects (i.e., enhanced spin 
temperature), or rapid motion in the obscuring gas. Ironically, the implications 
of these effects is that the obscuring disks must be small, typically not larger 
than a few tens of parsecs. 



Subject headings: galaxies: active: galaxies: Seyfert, galaxies; radio lines: 
galaxies 



1. Introduction 

The commonly accepted model explaining the power source for active galactic nuclei 
(AGNs) is a massive black hole fed by an accretion disk. The advantage of this model is the 
high efficiency of conversion of gravitational energy to radiant energy, but the main difficulty 
is the removal of angular momentum. Unlike stars, which constitute a dissipationless 
medium, gas and dust can lose orbital energy and angular momentum through collisions 
and shocks. At least in Seyfert galaxies, but probably in all AGNs, the host galaxy reserves 
a mass of neutral ISM comparable to or exceeding the lifetime needs of the AGN. Except 
in galaxy clusters, the mass of the hot phase component of the ISM is insufficient to fuel 
the most luminous AGNs over their duration. Although the transport mechanism remains 
a mystery, the conventional model holds that AGNs are powered by neutral ISM removed 
from the host galaxy (e.g., [Rees 19841) . 



Important variants on the standard model are the unifying schemes for AGNs. The 
dust associated with the neutral ISM can obscure a range of sight-lines to the AGN, 
affecting its outward appearance. It is proposed that both broad-line (e.g., Seyfert 1) and 
narrow-line (Seyfert 2) AGNs are identical, but the broad-line region (BLR) of Seyfert 2 
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galaxies may be obscured by a dusty ring, viewed more nearly edge-on, surrounding the 
active nucleus (this model is reviewed in, e.g., [Lawrence 1987| ; [Antonucci 1993| ; [Ward 1996| ). 



There are a handful of hidden Seyfert 1 galaxies discovered by spectropolarimetry (e.g., 
Antonuccl fc Miller l985| ; |lVliller fc Goodrich 199^ ; [I'ran, Miller, fc Kay 1991^1 ), but it is 



not clear whether the model is generally applicable. The detection of relatively high X-ray 
obscuring columns towards Seyfert 2 nuclei warrants the question of the nature and scale 
of the obscuring medium (e.g., [Awaki et al. 199"l| [Mulchaey, Mushotzky, fc Weaver 1992| ; 



Mulchaey et al. 1993| ; (Smith fc Done 1996| ; [Turner et al. 1997|) . 



We have observed a sample of 13 radio-bright Seyfert galaxies in H I (21 cm) absorption 
to explore the distribution and kinematics of the neutral ISM in the vicinity of the AGN. 
There are two primary advantages of this technique. First, absorption lines are potentially 
more sensitive than emission lines, the former limited primarily by the brightness of the 
background source rather than the column density or emission measure of the gas. Second, 
one can use powerful aperture synthesis arrays to resolve the absorption line gas on scales 
as small as the background radio source (i.e., potentially down to sub-parsec scales). 
Conversely, although the absorption experiment may not resolve the background source, 
the structure of the background source may be known from independent continuum studies, 
and so the effective resolution of the absorption experiment is improved. A third advantage 
is that one knows for certain that the absorption line gas must lie in front of the background 
source, removing the ambiguity suffered by emission line studies. Of course, the primary 
limitation is that one can detect and resolve absorption line gas only to the extent of 
the background source. In addition, because the opacity is inversely proportional to the 
excitation temperature, absorption lines are weighted more heavily by colder clouds. 

We present this work as follows. Section ^ describes the sample selection, observations, 
and data reduction techniques. Section ^ presents the data and detection statistics. 
In Section ^ we discuss the properties of the H I absorbing gas, including quantitative 
treatments of the spin temperatures and H I column densities. We also compare the 
observed H I columns with hydrogen columns based on X-ray spectroscopy. Section ^ 
describes the distribution and kinematics of the H I absorption based on the application of 
rotating disk models. Sections ^ and |^ discuss the implications of these observations on the 
standard model for AGNs, and we summarize our main conclusions in Section M. 



2. Observations and Data Reduction 
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2.1. Sample Selection 

The selection criteria for this survey were spiral galaxies with evidence for nuclear 
activity and radio brightness S^, > 50 mJy as resolved by the VLA A- Array. In essence 
this survey is a radio fiux-limited sample comprising mainly Seyfert galaxy nuclei. Our 
original survey consisted of 24 active spirals, but owing to a technical problem with our 
first observation (VLA program AB 605 on 19 Jul 1991), we were unable to recover the Hi 
spectra for these sources. We received a reduced time allocation for a repeat observation 
on 22 Nov 1992. Time and LST constraints reduced our sample to a subset of 13 of the 
brightest sources (VLA program AB 658). 

2.2. Observations 

We observed the sample galaxies with the VLA operating in its A-configuration, i.e., 
the longest baseline configuration, with baselines ranging from ~ 3.2-170 kA at 21 cm. 
The spatial resolutions in the image plane are typically ~ l'.'5-2". A summary of the main 
observing parameters are provided in Table |^. Receivers were tuned to the appropriately 
redshifted wavelength of the Hi 21 cm line (rest frequency u = 1420.406 MHz). We used the 
optical definition for velocity and referenced the velocities to the heliocentric frame. Eleven 
sources were observed in 4-IF mode, with paired IPs tuned redward and blueward of the 
systemic velocity to expand the velocity range covered by the receivers. The bandwidths 
are 3.1 MHz divided into 32 channels. In velocity units, the channel widths are ~ 21 km 
s~^, and, accounting for channel overlaps between the IF pairs, the total bandwidth is 
~ 1050 km s~^. The two remaining sources (NGC 2992 and NGC 5506) were observed using 
lA mode covering a single bandpass of 6.25 MHz, again with channel widths of ~ 21 km 

Source observations interleaved with shorter scans of a point-source calibrator for phase 
referencing. Each source observation was accompanied by a scan of a bright fiux standard, 
either 3C 48 or 3C 286. We also used these sources for bandpass calibration to remove the 
residual channel-to-channel complex gain variations. 

2.3. Data Reduction 



We reduced the data using standard procedures in AIPS. The fiux scale was set by 
scans of the radio continuum standards 3C 48 and 3C 286, defined according to |Baars et al.| 
(1977| ). Absolute calibration of interferometric phase was determined by point source models 
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for the phase reference cahbrators. We then employed 2-3 iterations of self-cahbration on 
each of our targets, properly including all of the confusing sources brighter than 5 mJy 
within two primary beam widths. The final CLEAN-deconvolved models were subtracted 
from the {u, v) database. We generated spectral-line cubes for each target by Fourier 
transform. The cubes were inspected for residual continuum emission or bandpass gradient 
not subtracted by the CLEAN-component model. For stronger sources, we subtracted a 
linear fit to the spectral baseline, but for weaker sources it sufficed to remove a continuum 
mean. All residual subtractions were performed in the image (i.e., not u,v) plane. 

A possible source of error in this experiment is a bandpass calibration artifact resulting 
from temperature variations, and hence length variations, of the waveguide system at the 
VLA ( [Carilli 1991| ; [Lilie 199^ ). The signature is a sinusoidal ripple with a period of 3 MHz, 



matching the bandwidth used in this experiment. Among the sources observed for this 
experiment, only the nuclear spectrum of NGC 1068 displayed broad spectral features that 
could be explained by the 3 MHz ripple. To test for this artifact, we generated a spectrum of 
the phase calibrator for NGC 1068. The phase calibrator showed similar spectral structure, 
corresponding to a spectral dynamic range of ~ 400. We therefore have little confidence in 
the broad spectral features of the nuclear 21 cm absorption profile of NGC 1068. 



3. Results 

We detect Hi absorption towards 9 out of the 13 sources observed in this survey, 
making a detection rate of 69%. The breakdown by nuclear activity is 3 out of 4 Seyfert 
1 nuclei detected (Fig. |I|); 4 out of 6 Seyfert 2 nuclei detected (Fig. and 2 out of 3 
starburst /LINER nuclei detected (Fig. H). Explicitly, the non-detections are Mkn 668 
(Seyfert 1), Mkn 3 and Mkn 348 (both Seyfert 2 with hidden Seyfert 1 nuclei), and 
NGC 2639 (LINER). We summarize the basic properties of the integrated Hi absorption in 
Table |^. Parameters of the integrated absorption line profiles are provided in Table ^. 

In several sources, the radio continuum was sufficiently resolved that we could isolate or 
map the Hi absorption distribution against the extended continuum. Sources in which the 
Hi absorption is unresolved, but the radio continuum is resolved, are NGC 2110 (Fig. 
NGC 2992 (Fig. 0), NGC 4151 (Fig. 0), and Mkn 6 (Fig. [n|). The continuum and Hi 
absorption in NGC 1068 are resolved (Figs. ^and|^), and the Hi absorption in NGC 3079 is 
marginally resolved against the background continuum (Figs. ^ and P). Follow-up MERLIN 



observations have better resolved the Hi absorption towards NGC 3079 ( Pedlar et al.| 



T996D , NGC 4151 (|Mundell et al. 1995|) , and Mkn 6 ( |Galhmore et al. 1998|) . Hi absorption 



towards Mkn 231 has been resolved by the VLB A ( [Carilli, Wrobel, fc Ulvestad 1997|) . 
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Appendix ^ details more general characteristics of these sources. 

The only published survey comparable to the present one is Dickey's (1986) non- 
imaging VLA survey of radio bright spiral and irregular galaxies. Accepting the small 
number statistics, our detection rate (69%) is better than Dickey's (53%). The most 
likely explanation is a ~ 4x improvement in sensitivity rather than sample selection. For 
example, our survey contains a greater fraction of AGNs relative to star-forming nuclei, but 
isolating the subsample of AGNs within Dickey's survey does not increase the detection 
rate. Explicity, nineteen sources were surveyed in pickey (1986|) , including 11 Seyfert or 
LINER nuclei and one BL Lac; rapid star-formation may dominate the radio continuum 
emission in the remaining sources. Dickey detected 10 out of 19 (53%) galaxies, or < 78% 
accounting for 5 additional, but marginal, detections. The statistics are comparable for 
the subsample of Seyfert and LINER nuclei: 5 (+4 marginal) out of 12 were detected, 
for a detection rate of 42% (< 75%). Furthermore, the distribution of continuum source 
brightnesses are similar between the present survey and [Dickey (1986|) . The improved flux 
sensitivity of the present survey constitutes an improvement in optical depth sensitivity. 
An example is the starburst nucleus NGC 3504, which was detected marginally in Dickey 
1986|) but is clearly detected in our survey. 



4. Properties of the H I Absorption Line Gas 



4.1. Spin Temperatures 

The column density Nhi and spin (excitation) temperature Ts are related to the 
observed Hi absorption spectrum according to: 



HI 



0.182 X 10^^ cm"2 



(T5/IOO K) 



Tt, dv 



where v is in units km s~^, and r„ is the opacity measured against the background 
continuum. In order to derive Hi column densities, we need some handle on the excitation 
state of the absorbing gas. 

There are three factors controlling the spin temperature of neutral hydrogen (|Field| 



1959a; Field 1959b; Field 1959c): collisions, pumping by Ly-a photons, and the absorption 



of 21 cm continuum radiation. The spin temperature approaches the kinetic temperature 
Tk in collisionally dominated (high density) regions. Hi absorption columns through dense 
clouds in galaxy disks typically match the Hi emission columns for a spin temperature of 
Ts ~ 50-300 K, close to the expected Tk (e.g., pickey, 'lerzian, fc Salpeter 1978| ; [Payne 



Salpeter, fc Terzian 1983| ; pijiszt et al. 1983] ; Praun 1997[ ). However, radiative effects are 
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expected to be important in the neighborhood of an AGN, ( Bahcall fc Ekers 1969 ), and Tk 
is a lower limit for T5. 



Fig. demonstrates the effects of 21 cm absorption and Lj-a pumping of the ground 
state of Hi in the 100 pc scale environment of a Seyfert nucleus (compare with Fig. 5 
O'Dea, Baum, fc Gallimore 1994]). For illustration, we modeled the 21 cm radiative 



m 



effects using the nuclear properties of NGC 4151. The background radio sources have a 
brightness temperature ~ 10® K at 21 cm (e.g., [Harrison et al. 1986|; [Ulvestad et al. 



1998|) , and the specific luminosity of the continuum near Ly-a is Fiya ~ 1.8 x 10^'' ergs s ^ 



Hz~^ ( [Kriss et al. 1992] ). In computing the density of Ly-a photons we further assumed a 
distance of 50 pc between the continuum source and the absorbing Hi cloud. A distance 
of 50 pc usually falls within the range of the observed, deprojected distances between the 
radio nuclei and Hi absorbing regions (see §^). We find that, in the 100 pc vicinity of an 
AGN, Ly-a pumping or 21 cm absorption controls T5 for densities N^i ^ 1000 cm~'^, but 
Ts settles asymptotically to Tk for nni ^ 1000 cm~^ as collisions dominate the excitation 
of the ground state. 



Note that Fig. |T2| is representative for typical Seyfert jets. A few Seyferts, however, 
have a bright self-absorbed core; Mrk 348 is an example. The intense radiation field near a 
self-absorbed radio core strongly affects Tgpin (e.g., Bahcall & Ekers 1969), and this effect 
may suppress H I absorption from an obscuring torus. We consider this possibility in more 
detail in 36. 3|- 



For most of the Seyferts in this survey, the Hi absorbed regions cover radii comparable, 
in projection, to the size of the NLR, and Hi absorption might plausibly arise from 
neutral filaments embedded within or shadowed by NLR plasma. Since oc {Ts Av)~^, 
Hi absorption studies are particularly sensitive to gas that is cold both thermally and 
kinematically. We would therefore expect most of the Hi absorption to arise from the 
coldest regions of clouds where hydrogen is primarily atomic. In AGN-irradiated clouds, 
the kinetic temperature of the cold atomic region is ~ 100 K ( [Maloney, Hollenbach, fc 



Tielens 1996|) . If we assume pressure equilibrium and typical NLR conditions Tg ~ 15000 K 
and Tie ~ 2000 cm~^ ( |Koski 1978 ), the density of the cold atomic region is uhi ~ 10^ cm"'^. 



At this limit collisions dominate, and Ts Tk (Fig. 0). 

We conclude that Hi absorption from the NLR region of Seyfert galaxies probably 
arise from the cold dense regions where Ts ~ T^ ~ 100 K. This result is remarkable in that 
Ts in cold Galactic clouds is also ~ 100 K, and so the conventional assumption T5 = 100 K 
holds up under a broad range of conditions. Of course, there may be situations where H l 
absorption traces a warmer atomic medium in the AGN environment; for example, the 
sight-line to the radio source may not cross the cold neutral region. In the warm atomic 
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region (and T5) can elevate to several 1000 K before hydrogen is significantly ionized 
( [Maloney, HoUenbach, fc 'I'ielens 1996| ). The assumption of T5 = 100 K therefore provides 
a lower limit to the true H I column. 



4.2. Column Densities 

We derived the column densities of the absorbing H I according to Eqn. |1] and assuming 
a spin temperature T5 = 100 K. The measured columns are listed in Table |^. The 
absorption opacity r was measured against channel-averaged continuum images from the 
same data set. Confusion between continuum sources can dilute the measured opacities or, 
equivalently, the Hi absorbing gas might not completely cover an unresolved background 
source. The reported values of Nhi are therefore lower limits to the true absorbing column. 
The columns for Mkn 6, NGC 4151, and Mkn 231 are derived from published MERLIN and 
VLBA observations ( |Mundell et al. 1995| ; |Carilh, Wrobel, fc Ulvestad 1997| ; [Gallimore eTaT 



19981) . 



One concern of this study is whether the H I absorbing gas is peculiar to the small-scale 
AGN environment or arises from the disk of the host galaxy on larger scales. Addressing 
this question, a simple prediction would be that the H I absorption should show no trend 
with inclination, since AGN disks appear to be oriented randomly with respect to outer 
galaxy disks (e.g., ^chmitt et al. 199^ ; |Nagar et al. 1999| ). Fig. |l^ compares the distribution 



of H I columns with the inclination of the host galaxy, estimated from optical axial ratios 
tabulated by |de Vaucouleurs et al. (1991| ) (compare with Fig. 5 in pickey 1986| ). Fig. |13 



also includes Hi absorption measurements from this work and Uickey (198(j| ). Only those 



sources with well-defined disk morphology and optical axial ratio are plotted. For reference, 
we also plot the range of Galactic Hi absorbing columns observed towards extragalactic 
radio sources ( [Payne, !Salpeter, fc 'lerzian 1982| ) and the radio source Sgr A in the Galactic 



Center ( pchwarz, Ekers, fc Goss 1982|) . 

We can draw two conclusions from Fig. |1^. First, running counter to the prediction 
for AGN disks, there is evidence for a trend in Hi column with host galaxy inclination. 
For instance, all four sources with i > 60° are detected in Hi absorption. We computed 
correlation statistics for the regression of N^i on seci using standard survival analysis 
techniques to account for the non-detections. The results indicate a real trend in nni with 
i: the probability for no correlation is 2.7% for Kendall's generalized r, < 1% for Cox's 
proportional hazard model, and 4.6% for Spearman's generalized p (although we caution 
that Spearman's p is not reliable for sample sizes ^ 30). This result suggests that in many 
Seyferts H I absorption arises in gaseous disks at least aligned with, if not cospatial with. 
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the kiloparsec-scale galaxy disks. We caution that the trend is only valid for the quantity 
Nni/Tspin, but a significant spread in T^pm would reduce confidence in the trend of Nhj 
with inclination. Nevertheless, we demonstrate in §^ that the kinematics and distribution 
of the H I absorption line gas are also consistent with the observed trend. 

The second conclusion is that the typical absorbing columns towards detected Seyfert 
galaxies exceed the Galactic columns by factors ranging from ~ 3-30, although there is 
some overlap for individual Seyfert galaxies and Galactic absorbing clouds. The mean 
face-on column for detected Seyferts is 4 x 10^^ cm~^, exceeding Galaxy columns by a factor 
of nearly 10. It would be interesting to compare the Seyfert Hi columns with a larger 
sample of non- active galaxies, but, lacking bright radio sources, non-active galaxies are not 
readily studied in Hi absorption. Nevertheless, based on Hi emission studies, the face-on 
columns through spiral galaxies are typically a few times 10^° cm~^ (e.g., surveys include 



the theses of |Bosma 1978| ; [Begeman 1987| ; [van Driel 198?1 ; their associated publications; and 



Proeils fc van Woerden 1994|) . Typical Hi emission columns of normal spiral galaxies are 



comparable to the absorption columns measured for the Galaxy rotated to a face-on view. 

The observed enhancement of H I column towards Seyfert nuclei is probably not an 
excitation effect. Demonstrated above in § [4.1| , the measured columns are lower limits 
because 21 cm absorption and Ly-a pumping in the vicinity of an AGN can raise Ts to 
1000 K or more. In contrast all of the Galactic values have been corrected for directly 
measured spin temperatures (for the Galactic Center column we assumed Ts ~ ~ 300 K 
in the 10 pc diameter disk after [Genzel et al. 1985| and references therein). Furthermore, 
owing to proximity, H I absorbing clouds in the Galaxy are better resolved compared to H I 
absorption towards Seyfert nuclei. Therefore, the Galaxy columns suffer less continuum 
dilution owing to confusion or covering fraction arguments. We conservatively conclude that 
the H I absorption columns towards Seyfert nuclei exceed the H I columns in the Galaxy, and 
perhaps other non-active spiral galaxies, by a factor of at least six (as shown in Fig. 



We show in §^ that H I absorption in Seyfert galaxies usually traces H I disks in the 
inner kiloparsec. In contrast. Hi columns are often reduced below detectability in the inner 
kiloparsec of nearby, non-active spirals (e.g., Roberts 1975| ; Giovanelli fc Haynes 1988 ). 



Limited by instrumental sensitivity, the central columns are typically <^ Cf 10^° cm~^, 
where C/ is the beam covering fraction of the Hi gas. The high detection rate for Seyfert 
nuclei argues for Cf approaching unity. Based on the interpretation of it appears that 
the H I surface density in the inner kiloparsec of Seyfert and starburst galaxies exceeds by 
a factor ^ 20 the H I columns in the centers of non-active galaxies. We discuss this result 
further in §|^. 
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4.3. Comparison with X-ray Columns 

X-rays provide another measure of the column of material obscuring the central engine 
in active galaxies. Photoelectric absorption by foreground gas that is not fully ionized 
attenuates the soft X-ray continuum spectrum of the central engine. The simplest models 
for AGN X-ray spectra usually involve a power law continuum modified by photoelectric 
absorption through the Galaxy and photoelectric absorption intrinsic to the source. The 
absorption is commonly parameterized by a total hydrogen column adopting standard cross 
sections and abundances (e.g., [Morrison fc McCammon 1983|) . 

Table ^ lists the X-ray absorption columns available in the literature for sources in the 
combined survey. There are often several models for a given source, and predictably the 
absorption column becomes poorly constrained as models become more complex. In each 
case, we selected the simplest model that provides a reasonable fit and the best constraints 
on an absorption column. Selecting models in this way tended to err on the side of caution; 
that is, more sophisticated models often give larger measures for the X-ray photoelectric 
column. The conclusions of this section are not affected by the selection of the simplest 
X-ray spectral models. 



Fig. compares the 21 cm absorption columns with the X-ray absorption columns. 
Under the hypothesis that the radio and X-ray nuclei are effectively cospatial, the prediction 
is a correlation between the 21 cm and X-ray columns. In addition, one would also expect 
that the X-ray columns should tend to be slightly greater than the 21 cm columns, since 
21 cm absorption is not sensitive to the presence of hydrogen locked in molecules or partially 
ionized gas. 

We find no measurable correlation between 21 cm and X-ray absorption columns, and 
the columns match reasonably only for Mrk 231. For most (8 out of 11) of the sources, the 
X-ray column exceeds the 21 cm column by at least an order of magnitude. We refer to 
these sources as having a "photoelectric excess" of absorption. There are three explanations 
for photoelectric excesses. First, it may be that most of the hydrogen is in molecular 
or ionized form. Alternatively, it may be that the primary X-ray absorber is dust, and 
the obscuring material has a lower gas-to-dust ratio than assumed in the modeling. The 
former explanation seems most likely for the Seyfert 1 nuclei, in which it appears that 
the absorbing medium is ionized (e.g., [George et al. 1998| ; [I'eldmeier et al. 1999[ ). The 
second explanation is that the radio emission may not be cospatial with the nuclear X-ray 
source, and so X-ray absorption and 21 cm absorption may probe different regions of the 
foreground media. We consider this possibility further in § The third explanation might 
be spin temperature effects. Increasing Tgpin to ~ 1000 K would produce a better agreement 
between X-ray and 21 cm columns for Mrk 6, NGC 2110, NGC 4151, and NGC 5506, and 
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the remaining photoelectric excess sources require Tspj„ > 10 K to explain the discrepancy. 
This explanation fails, however, to explain the result that H I absorption tends to avoid the 
nucleus in those Seyferts where the radio jet is resolved below). It seems more likely 
that 21 cm absorption does not trace the same gas giving rise to the soft X-ray opacity. 

The sources NGC 2992 and NGC 3079 show a "photoelectric deficit," meaning that the 
21 cm absorption column exceeds the X-ray photoelectric column. Again, it may be that 
the abundances assumed for the X-ray modeling are incorrect. For example, in this case, 
the gas-to-dust ratio may be higher than was assumed for the X-ray spectral modeling. 
The host galaxies of both sources are edge-on, and both show evidence for extended soft 
X-ray emission (e.g., |Elvis et al. 1990] |Kead, Ponman, fc Strickland 1997|; [Colbert, et al. 



The spatial resolution of the soft X-ray spectra is very coarse, ~ 30', much larger 
than the VLA beam. An alternative explanation, then, might be that extended X-ray 
emission dilutes the soft X-ray columns towards the nucleus. The extended emission that is 
resolved by current X-ray satellites is only ~ 10%-20% of the total, insufficient to explain 
the photoelectric deficits. Dilution of the soft X-ray columns would require extended X-ray 
emission on sub-arcminute scales not yet resolved; future X-ray satellites will have sufficient 
resolution to test the hypothesis. 

We conclude that there is no correlation between X-ray photoelectric columns and 
21 cm columns. Although there are a number of possible explanations, the simplest 
explanation may be that the radio and X-ray sources are not cospatial, and so 21 cm 
absorption and soft X-ray absorption do not (usually) trace the same foreground gas. We 
consider the implications in more detail in the following discussion (particularly below). 



5. Disk Models for the Distribution and Kinematics of the Hi Absorption 

There are now a handful of Seyfert nuclei in which the H I absorption is resolved against 
the background radio source: Mkn 6 (this study; [Gallimore et al. 1998|) , Mkn 231 ( |Carill 



IWrobel, fc Ulvestad 1997] ), NGC 1068 (this study. Fig. |; [Gallimore et al. 1994|) , NGC 2110 



(this study. Fig. NGC 2992 (this study. Fig. §, NGC 4151 (this study. Fig. |10D; [IVlundell 
[et al. 1995[ ), & NGC 5929 ( [Cole et al. T998[) . Hi absorption is also resolved in the starburst 
nucleus of NGC 3079 (this study; [Pedlar et al. 1996[ ; [^atoh, Inoue, k Nakai 19^5) ). The 
remaining detections that we have imaged, but in which the Hi absorption is not well 
resolved against the background continuum, are NGC 3504 and NGC 5506. 

It is remarkable that the Hi absorption tends to avoid the radio nucleus. More 
typically, the absorption is displaced by ;^ 100 pc in projection away from the central radio 
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source (Table |^). The exceptions are NGC 3079 and NGC 5506, which are edge-on, and 
NGC 4151, in which the Hi absorbing gas hes within a few tens of parsecs of the central 
engine ( [Mundell et al. 19951 ). We have demonstrated that detections are more likely with 



increasing inclination of the host galaxy (§[4 .21; Fig. O; see also pickey 1986|). Moreover 



the inferred absorption columns are not correlated with estimates of the X-ray absorption 
column (§ |4.3|) . It appears that Hi absorption in Seyfert and starburst nuclei traces gas in 
the inner disk of the host galaxy rather than gas associated with the AGN. In particular, 
the Hi absorption does not appear to be associated with a parsec-scale obscuring torus, 
with NGC 4151 being an important exception. 

To test this conclusion, we modeled the kinematics of the H I absorbing gas as disks in 
circular rotation. We spatially resolve the kinematics only in NGC 1068 ( [Gallimore et al.| 
19941), and so had to fit disk models to global absorption line profiles rather than spatially 



resolved channel maps. The results are listed in Table E and displayed in Figs. [T5|-23 



5.1. Procedure 

The disk models were computed and fit as follows. First, we generated a background 
continuum model from the continuum CLEAN components. For NGC 4151 we constructed 
a continuum model based on VLBI observations ( [Ulvestad et al. 1998|) , which recover all of 
the continuum flux in the MERLIN maps around the region of the Hi absorption ([Mundell 



et al. 19951) . We next generated a model disk described by vertical column density Nh] 



rotational velocity Vrou velocity dispersion a^, fixed at cr„ = 10 km s ^, as is commonly 
observed in spiral galaxy disks (e.g., [Shostak fc van der Kruit 1984|); scale height h, fixed at 



10 pc, self-consistent with the measured rotation velocities and disk radii; and, particularly 
for edge-on galaxies, inner and outer annular radii r^n and Tout- The scale height h for 
NGC 4151 was fixed at 0.1 pc because of the small radial scale of the absorbing disk. We 
also fitted the rotation curve as a power law v oc r", but, in every case, a ^ 0.01, meaning 
an essentially fiat rotation curve over the measured region. 

Since the hypothesis to be tested is whether the H I absorption traces normal disk gas, 
the inclination i and kinematic major axis (p were fixed to those values observed in the 
outer disk of the host galaxy in question. For NGC 1068, we took into account the rotation 
of the kinematic axis by 45° between the outer and inner galaxy. The apparent cause of 
this rotation is streaming in the central stellar bar ([Gallimore et al. 1994|) . Lacking spatial 



information, except for NGC 1068, we further assumed that the absorbing Hi uniformly fills 
the disk. In each case where the H I absorption is displaced from the nucleus, we assumed 
that the near side of the disk lies on the same side of the nucleus as the offset absorption. 
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The (deprojected) extent of the observed Hi absorption determines the inner and outer 
radii of the disk. For the edge-on galaxies in our survey, only the line profile constrains the 
inner and outer radii of the absorbing disk. 

Having generated and properly oriented the model disk and background radio source, 
we then integrated the line-of-sight opacity t{v) and computed an absorption spectrum for 
each background CLEAN continuum component. The final spectrum was created by adding 
up all of the model absorption spectra. We fit the model spectrum against the observed 
Hi absorption spectrum using a downhill simplex algorithm to minimize (chapter 10.4 
of [Press et al. 1992| and references therein). The best-fit model was then checked for 



self-consistency against the observed H I absorption morphology. We performed this check 
to ensure that the model disk would not obscure unabsorbed radio structure. 

After convergence, we estimated standard errors using a Monte Carlo technique. We 
assumed that the best-fit model line profile was the correct profile, and added random, 
Gaussian noise to the model spectrum to simulate an observed spectrum. The r.m.s. of 
the simulated noise was set equal to that of the residual spectrum (i.e., data — model). 
We repeated the fitting procedure on a series of ten simulated line profiles. The standard 
errors were derived from the deviations of parameters fitted to the simulated spectra. The 
failing of this technique is that it does not properly account for real absorption not present 
in the model and so may underestimate the real uncertainties where there are significant 
residuals. On the other hand, the measured standard errors appear to reflect the quality of 
the models and may not be far from their true values. 

The models are not unique in the sense that one can always contrive a model where all 
of the gas is falling into or flowing away from the background radio source. On the other 
hand, the structure of the background radio source leaves a "fingerprint" on the absorption 
profile that promotes the uniqueness of the fitted parameters of the disk model. That is, the 
structure of the background source will preferentially enhance or suppress the contribution 
of certain velocities to the spatially integrated line profile. An important failing is that 
our models cannot fit for spatially varying opacity in the absorbing gas because plausible 
opacity variations, such as spiral arm structure, require too many parameters for the fits to 
be meaningful. Nevertheless, the uniform disk model produces a reasonable and illustrative 
fit for all of the observed line profiles. 

One potential artifact of the fitting is a parsec-scale inner radius for the absorbing disk. 
This artifact results mainly from the need to fit the high velocity wings of the absorption 
line profile, since, for a given annular section of the disk at radius r surrounding a radio 
source of size tq, the velocity width is Av ~ Vrot[ro/\J (r — tq)^ + Tq]. The need for a small 
inner radius might be relaxed by allowing larger velocity dispersions or streaming velocities. 
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We chose to keep the dispersions fixed, since often the dispersions would grow unreahstically 
large in order to fit a discrete velocity component (e.g., NGC 1068, Fig. [16]; and NGC 5506, 
Fig. psp. Fitted values for the disk inner radius should be viewed as a lower limit. 



5.2. Results 

With the exception of NGC 4151, which is described further below, the disk models 
generally describe the observed Hi absorption profiles and locations very well. Poor values 
of can be traced generally to a discrete velocity component for which the simple disk 
model cannot account. Any discrete, residual components always arise within the velocity 
range traced by the model disk profile. We therefore suspect that unresolved opacity 
variations in the foreground disk may give rise to the discrete features in the line profile. For 
example, the spectrum of NGC 5506 shows a discrete redshifted component (Fig. |2^) that 
might arise from a cold cloud in the main disk which happens to fall along the line-of-sight. 

To demonstrate the validity of the galaxy disk model, we point to the peculiar 
absorption line profile of NGC 2110 (Fig. |T^). The absorption line is redshifted from 
systemic by nearly 300 km s~^, and the traditional interpretation would be infall towards 
the galaxy nucleus. However, the VLA images reveal that the absorbing gas is displaced 
from the nucleus along the kinematic major axis. Because of the special, chance alignment 
of the kinematic major axis and the radio jet, the northern lobe of the jet back-illuminates 
receding disk gas, and the result is a redshifted line profile. The uniform disk model 
reasonably reproduces the line profile shape, lending credence to this interpretation. The 
lack of a blue-shifted absorption line against the opposing, southern radio lobe indicates 
that the jet points in front of the disk to the south. 

NGC 4151 represents a special case, since the Hi absorption probably traces gas within 
the inner few tens of parsecs of the nucleus. Unfortunately, the location of the nucleus in 
the radio reference frame is not known, since all of the compact continuum sources are steep 
spectrum (e.g., comparing Pedlar et al. 1993| and Parral, Turner, fc Ho 1990|) . Since the 



Hi absorbed component (C4) is also free-free absorbed, and since the radio jet and NLR 
probably point toward us west of C4, the nucleus is probably near and west of C4 ( [Pedlar 



|et al. 19T?5| ). We selected the central VLBI component C4W as the location of the nucleus 
after [Mundell et al. (1995|). Discussed below, this placement only affects the quantitative 



results of the model but does not alter the qualitative conclusions. 

The radio jet of NGC 4151 does not align with the minor axis of the host galaxy, and 
so the AGN disk and galaxy disks are also misaligned. We therefore tried two disk models 
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for NGC 4151, one with an orientation matching the host galaxy ( Pedlar et al. 1993 ), and 
one matching the inferred properties of an AGN disk ( [Mundell et al. 1995| ). The AGN disk 
model fits the absorption spectrum better than the galaxy disk model. The reason is that 
the Hi absorption line centers very near the systemic velocity, and the AGN disk model 
places the H I absorbing gas along the disk minor axis. The galaxy disk model, on the other 
hand, predicts a slightly redshifted absorption component, since the absorption occurs 
towards the receding axis. As a result, the galaxy disk model requires counter- rotation to 
fit the observed Hi profiles. These conclusions hold so long as the radio jet originates at the 
kinematical center of the galaxy; the kinematical center may be placed anywhere along the 
radio jet and an AGN disk will always fit the absorption spectrum better than the galaxy 
disk model. 



5.3. Radial Streaming 

To look for evidence of infall that might fuel the AGN, we repeated the disk model 
fits allowing for a radial component to the gas velocity. The best-fit streaming velocities 
are listed in Table ^ The trivial, but unrealistic, solution would have been to set the 
rotational velocity to zero and assume purely radial motions. To avoid the trivial solution, 
we initialized the iterative fitting procedure with the parameters from the purely rotational 
model and allowed the circular and radial velocities to relax into a local minimum for x^. 

Although the error-bars are usually large relative to the fitted values, all of the 
streaming velocities are radially inward towards the nucleus except for NGC 1068 and 
NGC 4151. The deprojected infall velocities are typically a few tens of km s~^. 

Looking at the exceptions, the Hi absorbed region of NGC 1068 traces gas mainly on 
the trailing side of the central stellar bar ([Gallimore et al. 1994|) , over which the mean 
radial velocity of gas is expected to be away from the nucleus (e.g., [Roberts, Huntley, fc 



van Albada 1979| ; [Sanders fc Tubbs 1980|) . Although the streaming velocities in NGC 1068 



are very uncertain, the implied outward streaming is consistent with the predicted bar 
streaming. It is all but impossible to interpret the outflowing gas motions in NGC 4151 
since the absorption is not yet resolved against the VLBI structure. 

An important caveat is the main failing of these disk models; namely, we cannot 
recreate the distribution of absorbing g function of galactic coordinates. A localized 

enhancement of foreground column covering a finite range of radial velocities can mimic 
enhanced velocity dispersion or streaming velocities. Three-dimensional {a, 6, and z) 
models of the spectral line cubes would add enough information to allow at least a crude 
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treatment on the spatial distribution of H I absorbing gas, but the present data generally 
lack the spatial resolution for three-dimensional fits to be effective. The fitted streaming 
velocities can be viewed as upper limits to the true radial motions on ~ 100 pc scales. 

Accepting this caveat, the upper limits on the mass deposition rates constrain the 
luminosity that can be powered by the Hi disks. One concern, for example, is whether 
the upper limits on the infall rates might be insufficient to power the central engine. We 
have included in Table |^ the mass infall rates estimated by M ~ Mhi Vrad/fouu typical 
values are a few xO.l Mq yr~^. Among the Hi detections presented here, the mean, infall 
rate is 0.18 Mq yr~^. This rate can support a nuclear luminosity of L ~ 10^^ (e/0.1) ergs 
s~^, where e is the efficiency factor for mass-energy conversion in the accretion disk. This 
luminosity defines the traditional Seyfert/QSO luminosity boundary, and so the inferred 
mass deposition rates are comparable to that required to fuel Seyfert nuclei. We caution 
that this result does not constitute proof of infalling gas on 100 pc scales. On the other 
hand, we may conclude that the kinematics of the H I absorption are compatible with the 
standard model for AGNs, namely, that the active nucleus is fueled by gas removed from 
the host galaxy disk. 

We should point out that the nature of the H I absorption in active spirals appears to 
be very different from the Hi absorption in active ellipticals. Apart from a substantially 
lower detection rate, [van Gorkom, et al. (1989|) found that the Hi absorption lines were 



never blue-shifted, indicating either tangential motion or radial infall towards the central 
radio source. In contrast, we find that the kinematics of Hi absorption in Seyfert galaxies 
appear to be dominated by rotation, with a weak infall component superimposed. 



5.4. Stability Consideration 

Rotating gaseous disks are unstable to internal collapse and star formation where the 
disk exceeds a critical surface density S^. Neglecting, for now, the contribution of molecular 
gas and ionized gas to the local surface density, the ratio of the neutral hydrogen surface 
density T^hi to the critical density is 

~ Qm%N2lVrot/rpc , (2) 

where A''2i is the neutral hydrogen column density in units 10^^ cm~^, Vrot is the rotational 
velocity in km s~^, and is the radius in parsecs (based on P?oomre 1964| ; [Cowie 1981| ; 



[Kennicutt 1989|) . We measure T^hi/^c ~ 0.01-0.1 for the model Hi absorbing disks, based 
on the columns reported in Table The extreme values are <^ 0.001 for NGC 4151 (AGN 
disk) and NGC 5506 and 0.45 for NGC 4151 (Galaxy Disk). 
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The H I surface densities generally fall well below the critical density, but the surface 
density of molecular and ionized gas might still push the disks to instability. In the limit 
where atomic gas is the minority phase, the stability criterion becomes 

^\^] ■ (3) 



That is, for these disks to be unstable to star formation, the molecular and ionized gas 
masses would have to be 10-100 times the atomic gas mass over the inner kiloparsec. 

In global averages, the surface densities of atomic and molecular gas in spiral galaxies 
are similar to within factors of a few (e.g., [Kennicutt 1990| and references therein). On the 



other hand, gas in the centers of spiral galaxies tends to be predominantly molecular (e.g., 
Jackson et al. 1996| ; |Helfer fc Blitz 1997] ). Taking NGC 1068 as an example, [Heifer fc Blitz 



(1995| ) estimate a molecular gas column '^h2 ~ 5 x 10^^ cm ^ near the region we detect 



H I absorption. The molecular surface density exceeds the H I surface density by a factor of 
^ 20, pushing S/Ec 'i:^ 1. On the other hand, the rotation curve remains poorly constrained 
on this scale; a radial gradient of the rotation velocity would raise Sc and stabilize the disk. 

In NGC 1068, at least. Hi is the probably the minority phase within the inner 
100 pc. If other Seyfert galaxies are similarly rich in molecular gas, the disks traced by H I 
absorption disks may be unstable to star formation. Unfortunately, we are unable to test 
this suggestion with the present data. Current and future millimeter surveys, such as that 
currently taking place at the BIMA telescope (e.g., [Heifer et al. 19981) , may better address 
the stability issue in other Seyfert galaxies. 



6. Implications for Unifying Schemes 

Perhaps the most surprising result of this survey is the lack of H I absorption towards 
the AGN in Seyfert nuclei, particularly those known to harbor a heavily obscured central 
engine. The trivial explanation is that the central X-ray sources are radio-silent. However, 
the well-known hidden Seyfert 1 nuclei are also the most luminous Seyfert galaxies at 
radio wavelengths ( Moran et al. 1992 ), and it seems unlikely that the most luminous 



radio Seyferts contain radio-silent AGNs. Another possibility is that the obscuring disks 
are entirely molecular, or, alternatively, our sight-line crosses mainly through an ionized 
region of the disk. The line-of-sight X-ray columns commonly exceed 10^^ cm~^ (Table ^). 
Regardless of the viewing geometry, at some distance from the central ionizing source 
the obscuring material will become largely neutral. Photodissociation models predict a 
substantial atomic column detectable in 21 cm absorption, even if most of the material is 



- 18 - 



molecular; for example, in the presence of an Seyfert nucleus, dissociative excitations can 
liberate atomic gas within molecular clouds located out to distances of hundreds of parsecs 
( [Maloney, HoUenbach, fc 'I'ielens 19961 ). 

We next consider three ways in which H I absorption from AGN obscuring disks might 
be suppressed, namely kinematical line broadening, free-free absorption of near-nuclear 
radio continuum, and excitation effects for the 21 cm line. In each case the implication 
is that, ironically, the obscuring disk must be small, probably parsecs in scale, in order 
to suppress Hi absorption from the obscuring material. On the other hand, ^ 100 pc 
obscuring disks would not suffer these suppressing effects, and the prediction would have 
been strong H I absorption, contrary to the observations. 



6.1. Kinematical Line Broadening 

For a given mass of obscuring H I, the absorption line strength is inversely proportional 
to the velocity width of the line. In the environment of a compact mass, the absorption 
line-width should vary inversely with the (inner) radius from the AGN, and lines arising 
from rapidly moving gas near an AGN may be broadened below detectability. Without loss 
of generality, we consider a model where the ensemble of obscuring gas clouds is supported 
by random motion against gravitational collapse. Our 21 cm absorption experiment would 
be sensitive to columns: 

N„, « 3 X 10- c,„- r-^^'H,r- (g) (g^) . (4) 

where rpc is the distance from the AGN, Mj is the mass of the black hole in units 10^ Mq, 
(Tmjy is the channel map sensitivity in units of mJy, and Cioo is the background continuum 
strength in units of 100 mJy (e.g., as is appropriate for Mrk 348). We have also assumed 
that the local spin temperature is enhanced to ~ 8000 K, as may obtain in the warm 
atomic layer near an AGN ( [Maloney, Hollenbach, fc Tielens 199^ ). We conclude that 



the absorption experiment remains sensitive to column densities that may be typical for 
obscured Seyfert 1 nuclei, Nh ^ 10^^ cm~^, in the presence of kinematical line broadening. 
The obscuring gas would have to be within ~ 1 C{qq milliparsec of the AGN to suppress H I 
absorption below detectability. 
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6.2. Free-free Absorption 

The hypothesis in this case is that, as has been verified for NGC 1068 (e.g., |Gallimom 



Baum, fc U'Dea 1996| ; [Galhmore, Baum, fc U'Dea 199'^ ; |Koy et al. 19981 ), the nuclear 



continuum source has been attenuated by free-free absorption. Therefore, Hi in the neutral 
region of the obscuring torus lacks a background source on the smallest scales, and Hi 
absorption would not be detected. In this scenario, nuclear radio continuum does not trace 
emission from the AGN proper but instead from compact sources in the central, unresolved 
region of the radio jet. This is certainly the case for NGC 1068 (e.g., [Gallimore et al. 
11996 bp and Mkn 231 QUarilli, Wrobel, fc Ulvestad 1997] ), for example. 



Further motivating this hypothesis, the X-ray photoionization models of [Krolik fc 



[Lepp (19"89| ) and [Neufeld, Maloney, fc Conger (1994 ) have demonstrated that X-ray heated 



plasma within a parsec-scale obscuring torus can be optically thick to free-free absorption 
down to ~ 1 cm wavelength. Observationally, VLBI surveys are revealing free-free absorbed 
radio sources in other Seyfert nuclei (e.g., [Ulvestad et al. 1999|) . Pedlar et al. (1998|) 



find evidence for low-frequency free-free absorption towards the nucleus of NGC 4151. 
The paucity of obvious radio cores in Seyfert nuclei (e.g., [Wilson 1991| ), even relatively 



radio-bright Seyfert nuclei, also argues for free-free suppression of the Seyfert radio nuclei. 

We next demonstrate that the condition of free-free absorption places limits on the 
scale of the free-free obscuring plasma, which in turn limits the inner radius for neutral gas 
in the obscuring torus. For a given free-free opacity (i.e., emission measure), the free-free 
radio continuum fiux and hydrogen recombination line fiux scale with the projected surface 
area of the free-free absorbing plasma. Looking first at radio continuum emission from the 
nucleus, the constraint on the nuclear free-free radio surface brightness gives: 

r ^ 18 pc D20 a'Jl T^'" C]'/' , (5) 

where D20 is the distance in units of 20 Mpc, o"mjy is the channel map sensitivity in mJy, T4 
is the electron temperature in units 10^ K, and Cf is the covering fraction of the obscuring 
medium. This limit also assumes that the brightness temperature is equal to the electron 
temperature, appropriate for optically thick thermal plasma. 

Since, at 1" resolution, any free-free emission from the obscuring torus would probably 
be confused with synchrotron emission from the jet, we are limited not by the continuum 
sensitivity but rather the Hi absorption line sensitivity. (As new VLBI maps of Seyfert 
galaxies are published, continuum constraints on putative free-free absorbing regions may 
tighten the constraints argued here). X-ray absorbing columns of 10^'^ cm~^ (e-g-. Table |^) 
produce an Hi opacity of t^j ~ 20A23/(Tioo Afioo), where Tioo = T5/IOO K and /S.viqq 
is the line- width in units of 100 km s~^. A nuclear free-free continuum source would 
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be completely attenuated over the 21 cm line profile, and the limit on the AGN radio 
continuum flux is therefore the detection limit for a single channel, or roughly 3 mJy for 
this experiment. Cj ^ 50% based on the ratio of Seyfert 1 to Seyfert 2 galaxies ( |Lawrenc^ 
1991| ). Applying these values, the obscuring plasma must be smaller than r ~ 44 D20 pc, or 



else the background radio continuum provided by thermal free-free emission would produce 
an observable Hi absorption line. 

We next consider the hydrogen recombination line luminosity of the free-free obscuring 
region. We choose to consider hydrogen recombination lines mainly for ease of computation; 
consideration of forbidden line emission would require a photoionization model, which 
is beyond the scope of this work. Without loss of generality, we consider only Ha line 
emission. Since, according to unifying schemes, the ionized region of the torus would not be 
obscured by dust in Seyfert 1 galaxies, the spatially-integrated Ha luminosities of the NLR 
provide an upper limit to the emission from the obscuring region. For Seyfert galaxies, the 
median NLR Ha luminosity is Lho ~ lO'*^ ergs s~^ (based on characteristic line ratios in 
Ferland fc Osterbrock 1986| and observed H/3 luminosities from [Whittle 1992|) . Taking the 



usual assumption of Case-B recombination, the size of the free-free obscuring region must 
be: 

r ^ 28 pc (LhJIO'') C]"' rj}'' T^'-''' , (6) 

where Tff is the free- free opacity at the 21 cm rest wavelength, and we have used the usual 
radio approximation for the free- free opacity law from [Mezger fc Henderson (1967|) . 



If free-free suppression of the nuclear continuum source commonly occurs in Seyfert 
nuclei, as is at least suggested by this Hi absorption experiment, the implication is that the 
obscuring torus cannot be much larger than a few tens of parsecs. Turning the argument 
around, we would not predict substantial free-free absorption from a 100 pc scale obscuring 
torus; rather, we would have expected strong Hi absorption against a bright continuum 
source. 



6.3. Excitation Effects 

The third way to suppress Hi absorption is to raise the spin temperature, dropping the 
opacity for a given foreground column. We have demonstrated that, for typical radio source 
luminosities, the 21 cm transition in Seyfert obscuring disks is likely to be thermalized, and 
Hi absorption should be detected (§ [4.1| ). On the other hand, particularly luminous nuclear 
radio sources will drive the local spin temperature to the equivalent brightness temperature 
of the local 21 cm continuum radiation fleld (e.g., [Bahcall fc Ekers 19691). The radiation 



fleld from a compact nuclear source will drop as r ^, and so excitation effects will affect 
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21 cm absorption only if the obscuring gas is located near the central engine. Since the 
hypothesis is that we see the bright radio nucleus directly, this case can be viewed as the 
opposite extreme of the free-free absorption scenario presented above. 

In particular, spin temperature effects may explain the lack of H I absorption towards 
Mrk 348. Mrk 348 has an unusually bright nuclear radio source, Sy ~ 100-500 mJy at 
21 cm. This source is variable on time-scales of years, marking it as a likely location of the 
AGN (e.g., [Nell:' fc de Bruyn 1983| ; |Unger et al. 198^ . Assuming the source is much smaller 



than 1 pc, the 21 cm continuum temperature as a function of radius is 



Tr ^ 2.4 X 10^° K r'^ . (7) 



Applying the formalism of [Bahcall fc Ekers (1969| ), the 21 cm radiation field will dominate 



the effects of particle collisions at radii: 

r ^ 5 pc S\il nl'" D,, T^,^ , (8) 

where S'loo is the continuum source brightness in units of 100 mJy, is the gas density 
in units 10^ cm~^, Dgo is the distance to the source in units of 60 Mpc, and Tioo is the 
kinetic temperature in units of 100 K. The normalizations are specific to Mrk 348, but 
the result can be scaled and applied to any compact radio source. We estimate that the 
spin temperature might be as high as ~ 10^ K at a radius of 1 pc, and, as a result, 21 cm 
absorption opacity would be suppressed by a factor of 10^ T{^^\ 



6.4. Conclusions 



Sanders et al. (1989) originally proposed that the obscuring medium in Seyfert galaxies 
might span scales much larger than the conventional torus model. More recently, |Malkan,| 
Gorjian, fc Tam (1998 ) reported that optical images of Seyfert 2 galaxies are more likely to 



show ~ 100 pc dust absorption features than optical images of Seyfert 1 galaxies. They 
argued that these galactic dust lanes are primarily responsible for hiding the central engine 
and BLR in Seyfert 2 nuclei. The characteristic scale of the obscuring medium in this 
galactic dust model is ;^ 100 pc. To some degree galaxy disks must obscure Seyfert nuclei, 
since Seyfert nuclei are rarely found in edge-on disk galaxies ( |Keel 19g(]| ; |Simcoe et al. 

rwi). 



To explain the surprising lack of H I absorption towards known, hidden Seyfert 1 nuclei 
requires, ironically, obscuring disks on 10 pc scales or less. Turning the arguments around, 
we would expect that kinematical line broadening, free-free absorption, and excitation 
effects would have little effect on 21 cm absorption from gas on 100 pc scales. Therefore, 
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the galactic dust model predicts pronounced Hi absorption towards radio-bright Seyfert 
nuclei. We find instead that the dust lanes observed by [Malkan, Gorjian, fc 'lam (1998| ) are 
probably the source of the off-nuclear H I absorption reported here. We conclude that the 
inner scale of the obscuring medium in hidden Seyfert 1 galaxies must be less than a few 
tens of parsecs, in better agreement with the classical parsec-scale torus scenario. 

It is all but certain that free-free suppression of the nuclear continuum source has 
precluded the detection of 21 cm absorption towards NGC 1068. It seems very likely that 
excitation effects have suppressed 21 cm absorption in Mrk 348. The case of Mrk 3 is more 
difficult to address. [Kukula et al. (1999| ) identify a compact, flat spectrum VLBI source 
as the nucleus. The source is weak, only 9 mJy at 18 cm, and so excitation effects will 
be important only on sub-parsec scales (scaling from Eqn ^, above). On the other hand, 
this source might be intrinsically steep-spectrum but free-free absorbed at 18 cm, in which 
case it need not be identified as the nucleus. Additional spectral information is needed to 
interpret the nuclear radio source and the non-detection of 21 cm absorption. Otherwise, we 
have demonstrated that for plausible conditions in Seyfert nuclei, either free-free absorption 
or excitation effects would suppress 21 cm absorption from parsec-scale obscuring disks. 
These effects might naturally explain the tendency for 21 cm absorption to avoid the nuclei 
of obscured Seyfert galaxies. 

7. Implications for AGN Fueling 



One of the long-standing problems of AGN physics is the nature of the fuel source for 
the central engine. Based, for example, on source counts, or the sizes of radio jets, the 
lifetime of an AGN is expected to exceed ~ 10^ yr (e.g., as reviewed by Peterson 1997|) . 
Supporting a typical Seyfert luminosity of 10^^ ergs s~^ requires a fueling rate of about 
0.02 (e/0.1) Mq yr~^. Integrated over its lifetime, the central engine needs a total mass of 
~ 2 X 10^ (e/0.1) Mq. This mass is comparable to the central Hi masses derived in this 
absorption experiment (Table |^). The required infall velocities are low over the inner few 
hundred parsecs, only a few km s~^, and the Hi kinematics are compatible with our limits 
on the streaming velocities {vrad ~ few tens of km s"^; §|0|). We may also conclude that 
the inferred H l mass of the central galactic disks of Seyfert galaxies is sufficient to power 
the central engine over its lifetime. 

Bringing the gas to the central engine poses a more serious problem. The main obstacle 
is removing sufficient angular momentum to bring gas down to the inner parsec or so, where 
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viscous dissipation takes over (e.g., as reviewed by [Piiinney 1994|) . Stellar bars provide 
gravitational torque and may promote cloud-cloud collisions, funneling gas down at least as 
far as an inner orbital resonance. [Shlosman, i^'rank, fc Begelman (1989| ) proposed that mass 



accumulation at resonances can trigger the formation of nested bars. In principle, bars 
might nest down to arbitrarily small scales, carrying gas down to correspondingly small 
volumes. 

A basic test of this scheme is to compare the frequency of bars among active and 
non-active spiral galaxies. Taking advantage of new near-infrared imaging arrays, recent 
surveys have found that, in general, spiral galaxies are frequently barred (~ 70%), but the 
incidences of bars among active and non-active spiral galaxies are comparable ( [IViulchaey 



fc Regan 1997|). Moreover, the presence or strength of a bar appears to have little or no 



impact on the strength of the nuclear activity ([Ho, i^'ilippenko, fc Sargent 19971). As mass 



accumulates at the center of a bar, stellar orbits in the bar are perturbed, leading to the 
destruction of the bar (e.g., [Friedli fc Benz 1993| ; [Norman, Sellwood, fc Hasan 19961) . This 



mechanism might explain the absence of bars in a fraction of Seyfert nuclei. Nevertheless, 
the results of the imaging surveys raise two questions for the bar fueling scenario. First, 
there is growing evidence for black holes in the centers of both active and non-active galaxies 
(e.g., [Kormendy fc Richstone 1995[ ; [Miyoshi et al. 1995[ ; pallimore et al. 1996a[ ; [Eckart 



[fc Genzel 1996 ). Why are the black holes in barred, non-active spirals dormant? Second, 
other than the presence of a bar, what parameter might distinguish the host galaxies of 
active and non-active spirals? 

Our study suggests that the gas mass of the central galactic disk may be an important 
criterion distinguishing active and non-active spirals. The high detection rate of Hi 
absorption argues that the central, 100 pc scale disks of Seyfert nuclei are gas-rich, with 
azimuthally averaged surface densities S/^/(Seyfert) ;^ 6 (C/[Hl]/0.7) Mq pc~^, where 
C/[Hl] is the areal covering fraction of Hi over the central galactic disk, normalized to the 
detection rate. Such high surface densities should be easily detected in 21 cm emission in 
normal spirals. In contrast, non-active spiral galaxies commonly show an Hi depression 
over the inner few kiloparsecs (e.g., [Roberts 1975| ; [van Woerden, van Driel, fc Schwartzj 



1983[ ; [Giovanelli fc Haynes 1988|) , with column densities dropping below detections limits of 



~ 1 M0 pc-2. 

It is difficult to address this issue further, because there has been no systematic study 
of central Hi depressions in spiral galaxies. The reason is that, except in the dozen or so 
nearest galaxies, current aperture synthesis telescopes are generally not sensitive to Hi 
emission at sub-kpc resolution. To illustrate, we attempted to generate a comparison sample 
based on published Hi maps. The selection criteria were (1) Hubble type earlier than Sb, 
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since most Seyfert galaxies reside in early type spiral hosts; and (2) linear resolution better 
than 1 kpc, in order to reduce confusion with the outer disk. We chose galaxies from the 
Revised Shapley-Ames catalog ( [Sandage fc 'I'ammann 19811 ) and referenced them against 



the Martin's (1998) bibliographic index of 21 cm maps of galaxies. This selection process 
picked out 18 candidates for the comparison sample, but all but two of the candidates were 
either edge-on, active (Seyfert, starburst, or LINER), or strongly interacting with a close 
companion. This result reflects an observational bias towards the unusual but precludes an 
meaningful comparison between active and non-active galaxies. 

Nevertheless, we emphasize that feeding the AGN may require not only a stellar bar, 
but also an excess of gas mass in the inner kiloparsec as a source of fuel. Numerical studies 
of gas dynamics find that stellar bars remove a steady fraction /bar of the local gas mass 
per rotation; 0.01 ^ /bar ^ 0.1 ( [Athanassoula 1992| ; |Wada fc Habe 1992| ; [1^'riedh & Benz 



1993 ; Heller fc Shlosman 1994| ; [Athanassoula 1994 ). The gas mass required to fuel an AGN 



^7 T ( ^ \ / ./bar \ 



through bar-feeding is 

M = 1.01 X 10^ Me L44 (^— j i^roo, (9) 

where L44 is the luminosity of the AGN in units 10^^ ergs s~^, and fiioo is the bar pattern 
speed in units 100 km s~^ kpc^^. Accepting the considerable uncertainties in the efficiencies, 
the observed Hi absorption masses in active spirals, which give a lower limit to the total 
gas mass over the central disk, are comparable to the numerical predictions for bar feeding. 
Otherwise having a stellar bar and a central, supermassive black hole, non-active spirals 
may simply lack a central fuel source to power an AGN or nuclear starburst. Testing this 
scenario will require new Hi emission observations of a sample of active and non-active 
spiral galaxies, reasonably matched in Hubble type and physical resolution (distance). The 
prediction is that, over the inner kiloparsec, the mean Hi surface density of Seyfert galaxies 
will exceed the H I surface density of non-active spirals. 



8. Conclusions 

Our main conclusions are as follows: 

1. Hi 21 cm absorption in Seyfert and starburst galaxies appear to trace gas in the 
central disks of the host galaxy. The radial scale of the central absorbing disks 
is hundreds of parsecs (i.e., less than a kiloparsec). High- velocity Hi absorption, 
Doppler-shifted away from the systemic velocity of the host galaxy, can usually be 
explained by normal rotational motion rather than radial streaming (e.g., NGC 2110). 
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2. The Hi kinematics are compatible with, but do not require, streaming velocities of 
several tens of km s~^. The corresponding mass inflow rates can be several Mq yr~^, 
generally exceeding the mass transfer rate required to fuel a Seyfert nucleus. 

3. NGC 4151 appears to be a special case ( |lVlundell et al. 1995| ). The kinematics 
and distribution of the Hi absorption are consistent with a compact, parsec-scale 
disk oriented perpendicular to the radio jet. In contrast, the Hi properties are not 
compatible with the larger scale galaxy disk orientation and rotation sense. 

4. Hi absorption in Seyfert galaxies generally avoids the nucleus but is more commonly 
seen against an extended radio component. The 21 cm Hi columns are not correlated 
with the hydrogen column derived from X-ray spectroscopy, even in known hidden 
Seyfert 1 nuclei. We propose that Seyfert nuclei are commonly free-free absorbed, 
and no 21 cm continuum from the AGN escapes to the neutral layer of the obscuring 
medium. Alternatively, bright nuclear continuum sources, such as in Mrk 348, can 
enhance the spin temperature and suppress 21 cm hne absorption. These suppression 
mechanisms constrain the size of the obscuring medium to smaller than a few tens of 
parsecs. 

5. We suggest that the central disks of active spiral galaxies may show higher Hi surface 
densities compared with non-active spirals. This result may resolve the observational 
challenges for the bar fueling model of AGNs. That is, to be active a spiral galaxy 
might require both a central stellar bar and an excess of neutral gas. Testing this 
prediction will require further studies of 21 cm emission from matched samples of 
active and non-active spirals. It is not clear from this study whether we are witnessing 
gas en route from the outer galaxy (i.e., the fuel line) or a local reservoir from which 
the nucleus derives its fuel (i.e., the fuel tank). 
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A.l. Detected Sources 



Mkn 6 is cL Seyfert 1.5 galaxy. Hi absorption was originally reported by [Heckman, Balick 



|fc Sullivan (1978D and [Baan fc Hascliick (1983| ), whose detections we confirm. We have 
reported a more detailed follow-up MERLIN observation QGallimore et al. 199S ). The 



nuclear radio source, marginally resolved by the present observation, collimates into a 1" 
(360 pc) radio jet ([Kukula et al. 199^ ). Hi absorption is detected against the northernmost 
radio component of the 1" nuclear jet. The offset absorption appears to arise from a 
hundred-pc scale arm or dust lane crossing just north of the AGN. The AGN itself may be 
partially absorbed by the dust lane, which may lend to its classification as a Seyfert 1.5. 
Otherwise, there are no radio candidates for the AGN. Recent ASCA spectroscopy suggests 
a high emission measure to the nucleus, and so the radio nucleus may be free-free absorbed 
( P^'eldmeier et al. 1999| ). 



Mkn 231 is a luminous and morphologically peculiar Seyfert 1 galaxy. Hi absorption was 
first reported by [Heckman, Balick, fc Sullivan (19"78| ) and ([van der Hulst, Hummel^^ 



Dickey 1982[) . Our VLA spectra confirm these earlier detections. The VLBA observations 



of [Carilh, Wrobel, fc Ulvestad (1997[) resolve the H I absorption against a 0'.'44 (~ 350 pc) 



radio continuum disk surrounding the AGN. The H I absorption displays a ~ 220 km s ^ 
velocity gradient along the major axis of the diffuse radio continuum disk, and [Uarilli 



[Wrobel, fc Ulvestad (199"^ ) argue that the H I absorption arises in the atomic component of 
a molecular disk mapped by millimeter-wave CO emission ( [Bryant fc Scoville 19961) . 



NGC 1068 is the prototype of Seyfert 2 nuclei but is also the first Seyfert found to have a 
Seyfert 1 spectrum in polarized light ( [Antonucci fc Miller 1985[ ). [Baan fc Haschick (1983| ) 



originally reported Hi absorption seen as a sharply defined trough in the Hi emission profile. 
Our observations resolve out H I emission, and the H I absorption is spatially resolved. The 
main results were originally presented in pallimore et al. (1994[ ). The bulk of the absorption 
traces the inner disk of the host galaxy, seen against the southwestern radio structure. 
There is kinematic evidence that the gas is responding to the kiloparsec scale stellar bar. 
Closer to the nucleus the absorption hues broaden, and there appear multiple velocity 
components. The opacity of the absorption line gas also appears to increase, but towards 
the nucleus the opacity seems to drop sharply. VLBI studies have failed to detect the 
nucleus (radio component SI) in 21 cm radio continuum ( [Roy et al. 1998[ ). The apparent 
drop in opacity is probably a resolution effect owing to the lack of a background radio 
source. There may be a broad absorption profile very near the nucleus, whose kinematics 
are consistent with the rotating H2O maser disk ( [Gallimore et al. 1996a] ), but bandpass 



calibration uncertainties preclude a confident detection in the present observations. 



NGC 2110 is an X-ray bright Seyfert 2 galaxy ( [McClintock et al. 1979[) . At sub-arcsecond 
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resolution, the "figure-eight" radio structure of our 21 cm continuum image resolves into 
a lazy-S-shaped radio jet ( [Wilson fc Ulvestad 1983| ). We detect redshifted Hi absorption 



offset to the south towards the weaker radio lobe. 

NGC 2992 is a narrow line X-ray galaxy ([Ward et al. 19801) . The radio structure comprises 
an unresolved nuclear source surrounded by figure-8 loops ( [Wehrle fc Morris 1988[ ). Hi 



absorption is detected offset from the radio nucleus against one of the extended lobes. The 
host galaxy is viewed nearly edge-on, and so it seems likely that the H I absorption traces 
gas in the foreground disk. 



NGC 3079 is classified as a LINER ([Heckman 1980[) with evidence for active star-formation 



e.g., [Condon fc Broderick 198"8[ ; [Irwin fc Seaquist 19901 ). The Hi absorption was discussed 



briefly by [Galhmore et al. (1994]) but in more detail by [Baan fc Irwin (1995[ ). MERLIN 
observations well resolve the rotational trend of the absorption line gas ([Pedlar et al. 1996| ), 
which matches the sense of the rotation traced by H I emission ( ^rwin fc Seaquist 1991|) . The 
bulk of the absorption probably arises from the edge-on but normally rotating disk, observed 
in silhouette against extended, diffuse radio emission associated with star-formation ( jPedlar 



et al. 1996[) . Based on new VLBI observations, [Satoh, Inoue, fc Nakai (1998[ ) report that 



the H I absorption viewed towards the central VLBI radio jet appears to be counter-rotating 
relative to the gas motions in the outer galaxy, but it is unclear whether the H I absorption 
on VLBI scales traces gas in rotation as opposed to, say, nuclear-driven outflow. The 
apparent complexity of the broad absorption line proflle may reflect the structure of the 
VLBI jet as much as the kinematics and distribution of the absorbing gas. 

NGC 3504 is a barred starburst nucleus originally included on Sersic's list of peculiar nuclei 
( [Sersic 1973[ ). The nuclear radio emission comprises a compact, unresolved core surrounded 
by a ~ 10" diffuse halo (e.g., [Saikia et al. 199"^. The halo probably traces synchrotron 



emission from a 12" diameter starburst ring, revealed in near-infrared color maps ( [Elmegreen[ 
et al. 1997|) . Pickey (1986[) reported a tentative detection of Hi absorption. We clearly 



detect an absorption line towards the unresolved radio nucleus of this galaxy. 

NGC 4-151 is a well-studied Seyfert 1.5 galaxy. Hi absorption was originally detected in this 
experiment and followed- up by MERLIN observations ( [Mundell et al. 1995[ ). The nuclear 
radio structure resolves into a ~ 5" (~ 320 pc) radio jet ( [Pedlar et al. 199^ and references 
therein). Hi absorption is detected only towards the radio component C4 ( pVlundell et al 



[1995[) , which itself resolves into a double source on VLBI baselines ( [Harrison et al. 1986[ ) 



The tight conflnement of the Hi absorption to the VLBI portion of the radio jet argues 
that the absorption arises from near-nuclear gas, perhaps a disk surrounding the central 
engine, at tens of parsecs from the central engine ([Mundell et al. 1995[) . 
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NGC 5506 is an X-ray bright Seyfert 2 galaxy (|Wilson et al. 1976t [Rubin 1978| ). Hi 
absorption was originally reported by |Tliuan fc Wadiak (1982| ). The radio emission is 
dominated by a compact radio source, unresolved at ~ O'/l resolution ( [Wehrle fc Morris 
1987 ). A low surface brightness arc extends ~ O'.'S north from the nucleus ( |Unger et al.| 
19861 [Wehrle fc Morris 1987|) , and faint emission extends ~ 30" from the nucleus ([Colbert[ 



[et al. 19%D . The host galaxy is viewed nearly edge-on, and it seems likely that most of the 
absorption arises in foreground disk gas. 



A. 2. Non-detections 



Mkn c? is a Seyfert 2 nucleus harboring a hidden BLR ([Miller fc Goodrich 19901) . The radio 
structure is linear, resolving into a 2" (~ 600 pc) jet ( [Kukula et al. 19931) . Surprisingly, we 



do not detect Hi absorption in this hidden Seyfert 1. On the other hand, there are no fiat 
spectrum components that clearly mark the location of the central engine ( [Kukula et al 



1993[) . It may be that the nuclear radio source is free-free suppressed or radio silent, and so 



there is no background radio source to illuminate the material obscuring the central engine. 

Mkn 348 (NGC 262) is a Seyfert 2 galaxy that is unusual in two respects. Firstly, it has a 
very large H I disk grown perhaps by a tidal encounter with a neighboring spiral ( [Heckman 



[et al. 198^ ; [Simkin et al. 19871) . Secondly, there is evidence for a hidden BLR in this source 



(e.g., [Miller fc Goodrich 199Q . Although the nucleus hosts a very bright subarcsecond 
radio triple (e.g., [Netf fc de Bruyn 1983[ ), we detect no Hi absorption in this experiment. 



Halkides, Ulvestad, fc Roy (1997] ) report evidence for free-free absorption below 6 cm 
wavelength, but it is unclear which components of the radio triple might be absorbed. (This 
result is consistent with the spectral decomposition of [Netf fc de Bruyn 1983| ). The lack of 



H I absorption from the medium obscuring the BLR might again be ascribed to free-free 
suppression of the nuclear radio continuum source at 21 cm wavelength. 

Mkn 668 (OQ +208) has the optical spectrum of a Seyfert 1 but, given its radio power, may 
be more appropriately classified as a broad line radio galaxy (e.g., [Stanghellini et al. 1997[ 



and references therein). Mkn 668 is also notable for being one of the nearest GHz-peaked 
spectrum radio sources ( [O'Dea, Baum, fc Stanghellini 1991[ ). VLBI observations resolve the 
compact radio source into a 7 mas (~ 5 pc) double ( [Stanghellini et al. 1997[ ). 



NGC 2639 is a LINER hosting an H2O megamaser (jWilson, Braatz, fc Henkel 199^ ). As may 
also be the case for the megamaser sources NGC 1068 and NGC 4258, free-free absorption 
may suppress 21 cm continuum within the maser disk ([Gallimore, Baum, fc O'Dea 199"^ ; 
Herrnstein et al. 1998[ ). This may explain the lack of Hi absorption corresponding to the 
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H2O maser. Alternatively, there may be very little atomic gas associated with the H2O 
masers, but that seems unlikely in the photodissociative environment of an AGN. The radio 
emission resolves into a 1'.'4 linear structure ( |Ulvestad fc Wilson 19891 ), but, mainly because 
no spectral index information is available, there is presently no radio candidate for the 
AGN. 
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Fig. 1. — Global (spatially-integrated) H I absorption spectra of Seyfert 1 nuclei. The spectra 
are continuum-subtracted and corrected back to zero rest velocity based on the observed 
recessional velocities in Table ^ The only non-detection among the Seyfert 1 nuclei was 
Mkn 668. 
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Fig. 2. — Global Hi absorption spectra of Seyfert 2 and narrow-line X-ray nuclei. The 
spectra are continuum-subtracted and corrected back to zero rest velocity based on the 
observed recessional velocities in Table ||. The non-detections are Mkn 3 and Mkn 348. 

Fig. 3. — Global Hi absorption spectra of starburst and LINER nuclei. The spectra 
are continuum-subtracted and corrected back to zero rest velocity based on the observed 
recessional velocities in Table |l|. The only non-detection is NGC 2639. 

Fig. 4. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 1068. The contour levels are 0.21, -0.21, -0.41, -0.60, & -0.80 Jy km 
s~^. The maximum grayscale level is clipped at 0.3 Jy. The star marks the location of the 
nucleus. The beam is indicated by the filled ellipse in the lower left corner. 

Fig. 5. — Channel-map overlays of Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 1068. The contour levels are -13.9, -9.2, -6.1, -4.1, & -2.7 (3a) 
mJy beam~^. The grayscale stretch is logarithmic, and the maximum grayscale level 
displayed is clipped at 1.0 Jy. The channel velocities, in km s"^, are provided in the upper 
right corner of each panel. 

Fig. 6. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 2110. The contour levels are 0.45, -0.45, -0.57, -0.71, & -0.90 Jy km 
s~^. The star marks the location of the radio nucleus. The beam is indicated by the filled 
ellipse in the lower left corner. 

Fig. 7. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 2992. The contour levels are 0.36, -0.36, -0.49, -0.66, & -0.90 Jy km 
s~^. The star marks the location of the radio nucleus. The beam is indicated by the filled 
ellipse in the lower left corner. 

Fig. 8. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 3079. The contour levels are 0.57, -0.57, -1.1, -2.2, -4.3 & -8.4 
Jy km s~^. The continuum is displayed with a logarithmic stretch. The star marks the peak 
of the radio continuum, which is presumably the nucleus. The beam is indicated by the filled 
ellipse in the lower left corner. 
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Fig. 9. — Channel-maps of the Hi absorption (contours) towards NGC 3079. The contour 
levels are —40.5, —19.1, —9.0, —4.2, —2, & —2 mJy beam~^. The continuum is displayed 
as contours on greyscale in the lower left panel. The star symbol marks the location of 
the radio continuum peak on the channel maps. The continuum contour levels are —0.30, 
0.30, 0.66, 1.5, 3.2, 7.0, 15.4, 34.0, & 74.7 mJy beam~^. The channel velocities, in km s~^, 
are provided in the upper right corner of each panel. A ~ 0'.'3 scale velocity gradient, with 
velocities increasing from northwest through southeast, is marginally resolved by our data 
(compare with Pedlar et al. 1996). 



Fig. 10. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of NGC 4151. The contour levels are 0.3, -0.3, -0.45, -0.70, & -l.lJy km s~^ 
The beam is indicated by the filled ellipse in the lower left corner. 



Fig. 11. — Overlay of integrated Hi absorption (contours) on radio continuum (inverse 
grayscale) of Mrk 6. The contour levels are 0.19, —0.19, —0.46, —1.1, & —2.7 Jy km s~^. 
The star marks the location of the optical nucleus (Clements 1981). The beam is indicated 
by the filled ellipse in the lower left corner. 



Fig. 12. — Plots of T5(21 cm)/Ti^ for varying kinetic temperature Tk as a function of neutral 
hydrogen density nn- The models plotted in the bottom panel include only the effects of 
21 cm absorption and collisional excitation, but the plots in the top panel include the effects 
of Ly-a pumping of the ground state of Hi. The three curves on each panel correspond 
to Tk = 3, 10, and 100 K. We modeled the radio and UV continuum properties after 
NGC 4151, assuming a distance of 50 pc between the AGN and absorbing cloud (further 
details are provided in the text). Irrespective of the presence of a UV pump, the 21 cm 
transition thermalizes for Uh ^ 1000 cm~^. Note that 21 cm continuum excitation may 
affect Tgpin in Mrk 348 and other sources with bright, flat-spectrum cores (see § |6.3|) . 
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Fig. 13. — Absorbing Hi column density plotted vs. galaxy inclination. The columns 
for Seyfert galaxies are plotted as filled circles; upper limits are indicated by arrows. 
Measurements from Dickey (1986) have been included. The Seyfert columns assume a 
characteristic Tspin = 100 K, which may be a lower limit to the true spin temperature. The 
dashed line traces the secant curve appropriate for a face-on column density of 2 x 10^^ cm~^, 
representative of the Seyfert detections. The open circles mark Galactic absorption data 
from Payne, Salpeter, & Terzian (1982) and Schwarz, Ekers, & Goss (1982). The Galaxy 
absorption columns have been corrected for measured T^pin- The effective inclination of the 
Galaxy absorption measurements was taken to be the complement of the Galactic latitude. 
The mean face-on column for the Galactic absorption measurements is 3 x 10^° cm~^ and is 
traced by the solid secant curve. Compare with Figure 5 of Dickey (1986). 



Fig. 14. — A plot comparing Hi columns measured by 21 cm absorption (x-axis) with 
total hydrogen columns measured by model fits to soft X-ray spectra (y-axis). We assumed 
Tspin = 100 K in deriving Nhi] the Nhi measurements are actually lower limits. Narrow-line 
objects, including Seyfert type 2, LINER, and starburst nuclei, are marked by filled squares, 
and broad-line objects are marked by unfilled squares. The line traces Nhi{21 cm) — Nh{X.- 
ray). 
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Fig. 15. — A kinematic model fit for the H I absorbing disk of Mkn 6. The data are MERLIN 
observations from Galhmore et al. (1998). Left: A schematic of the Hi absorbing disk modeL 
The disk has been shaded so that the far side appears hghter than the near side, as if it 
were viewed through an atmosphere; the effect is to mimic the appearance of a dust lane. 
The disk shading of the is lightened where there is significant background continuum, mainly 
to improve the visibility of the continuum contours. We erased the disk shading where the 
radio continuum structure lies foreground to the disk to illustrate the relative orientation of 
the disk and jet. The scale height of the disk has been exaggerated for presentation. The 
contours trace the 20 cm radio continuum image. Note that the outer radius of the H I disk 
is a lower limit set by the extent of the observed Hi absorption. Right: Model spectra. 
The top panel shows the fit (thick line) against the Hi spectrum (thin line), and the bottom 
panel shows the residuals. 



Fig. 16. — A kinematic model fit for the Hi absorbing disk of NGC 1068. Left: A schematic 
of the H I absorbing disk model, following the convention of Fig. [1^. The contours trace the 
20 cm radio continuum image. The outer radius of the Hi disk is a lower limit set by the 
extent of the observed H I absorption. Right: Model spectra. The top panel shows the fit 
(thick line) against the Hi spectrum (thin line), and the bottom panel shows the residuals. 



Fig. 17. — A kinematic model fit for the Hi absorbing disk of NGC 2110. Left: Schematic 
view, following the convention of Fig. |T5[ The contours trace a 6 cm radio continuum image 
from Ulvestad & Wilson (1984). We display the higher resolution VLA image to emphasize 
the H I absorbed radio structure. The outer radius of the H I disk is a lower limit set by the 
extent of the observed Hi absorption. Right: Model spectra. The top panel shows the fit 
(thick line) against the Hi spectrum (thin line), and the bottom panel shows the residuals. 



Fig. 18. — A kinematic model fit for the Hi absorbing disk of NGC 2992. Left: Schematic 
view, following the convention of Fig. |15[ The contours trace a 6 cm radio continuum image 
from Ulvestad & Wilson (1984). The outer radius of the Hi disk is a lower limit set by the 
extent of the observed H l absorption. Right: Model spectra. The top panel shows the fit 
(thick line) against the Hi spectrum (thin line), and the bottom panel shows the residuals. 



Fig. 19. — A kinematic model fit for the Hi absorbing disk of NGC 3079. Left: Schematic 
view, following the convention of Fig. |T5[ The contours are 20 cm radio continuum. Right: 
Model spectra. The top panel shows the fit (thick line) against the Hi spectrum (thin line), 
and the bottom panel shows the residuals. 
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Fig. 20. — A kinematic model fit for tfie Hi absorbing disk of NGC 3504. Left: Schematic 
view, following the convention of Fig. 15. The outer radius of the disk is poorly constrained, 
because the background continuum is so weak, and therefore offset H I absorption contributes 
little to the global profile. The contours are 20 cm radio continuum. Right: Model spectra. 
The top panel shows the fit (thick line) against the Hi spectrum (thin line), and the bottom 
panel shows the residuals. 



Fig. 21. — A kinematic model fit for the Hi absorbing disk of NGC 4151. Left: Schematic 
view, following the convention of Fig. |T3|. The contours are 20 cm radio continuum modeled 
after the VLBI observations of Ulvestad et al. (1998). In this model, the disk is oriented 
with the same inclination and position angle as the outer galaxy disk. Right: The top panel 
shows the fit (thick line) against the Hi spectrum (thin line), and the bottom panel shows 
the residuals. Note that this model requires a counter-rotating disk. 



Fig. 22. — A kinematic model fit for the Hi absorbing disk of NGC 4151. Left: Schematic 
view, following the convention of Fig. |15|. The contours are 20 cm radio continuum modeled 
after the VLBI observations of Ulvestad et al. (1998). In this model, the disk is oriented 
with the same inclination and position angle as that inferred for an AGN disk (Mundell et 
al. 1995). Right: The top panel shows the fit (thick line) against the Hi spectrum (thin 
line), and the bottom panel shows the residuals. 



Fig. 23. — A kinematic model fit for the Hi absorbing disk of NGC 5506. Left: Schematic 
view, following the convention of Fig. |15[ The contours are 20 cm radio continuum. Right: 
Model spectra. The top panel shows the fit (thick line) against the Hi spectrum (thin line), 
and the bottom panel shows the residuals. 



Table 1. Source Properties 



Source 


Right Ascension 


Declination 


cz 




Ref. 


Activity 


Morpholog 






(1950) 




(1950) 




km s^ 


-1 






(NED)^ 


Mkn 348 


00 


46 


04.85 


+31 


41 


04.5 


4507 


±4 


3 


2(1) 


SA(s)0/a 


NGC 1068 


02 


40 


07.07 


-00 


13 


31.2 


1150 


5 


1 


2(1) 


SA(rs)b 


NGC 2110 


05 


49 


46.38 


-07 


28 


00.9 


2335 


20 


2 


2 


SABO- 


Mkn 3 


06 


09 


48.24 


+71 


03 


10.5 


3998 


6 


3 


2(1) 


SO 


Mkn 6 


06 


45 


43.96 


+74 


29 


09.9 


5640 


10 


7 


1.5 


SABO+ 


NGC 2639 


08 


40 


03.12 


+50 


23 


10.2 


3336 


11 


3 


L 


SA(r)a 


NGC 2992 


09 


43 


17.62 


-14 


05 


42.9 


2314 


6 


3 


X 


Sa pec 


NGC 3079 


09 


58 


35.01 


+55 


55 


15.5 


1124 


10 


4 


L 


SB(s)c 


NGC 3504 


11 


00 


28.53 


+28 


14 


31.3 


1534 


2 


5 


SB 


SAB(s)ab 


NGC 4151 


12 


08 


01.01 


+39 


41 


02.1 


997 


3 


6 


1.5 


SAB(rs)ab 


Mkn 231 


12 


54 


05.03 


+57 


08 


38.2 


12650 


20 


8 


1 


SA(rs)c? 


Mkn 668 


14 


04 


45.61 


+28 


41 


29.3 


22957 


40 


9 


1 




NGC 5506 


14 


10 


39.14 


-02 


58 


27.1 


1815 


9 


3 


2 


Sa pec 



Note. — Properties of the sources observed in this experiment. Coordinates mark the 21 cm 
continuum peak, not necessarily the nucleus, at ~ 2" resolution. References: (1) [Brinks et al 



(T997|) ; (2) [Nelson fc Whittle (1995| ); (3) [de Vaucouleurs et al. (199l[ ); (4) [Irwin fc Seaquist 



T99l[) ; (5) IKenney, Carlstrom, fc Young (1993[ ); (6) [Pedlar et al. (1993| ); (7) [MeaburnJ 
Whitehead, fc Pedlar (19891) ; (8) [Rigopoulou et al. (199q ); (9) [Huchra et al. (199q) . Activity 
codes: l.n = Seyfert l.ra; 2 = Seyfert 2; 2(1) = hidden Seyfert 1; L = LINER; SB = starburst; 
X = narrow-line X-ray galaxy (see notes in 

'^This research has made use of the NASA/IPAC Extragalactic Database (NED) which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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Tablc 2. Observing Parameters 



Source 


tint 


RMS 


Continuum 




Beam 






(min) 


(mJy beam^^) 


(mJy beam~^) 




(arcsec) 




Mkn 3 


20 


1.5 


895 


2.3 


X 1.6 


72° 


Mkn 6 


52 


1.0 


176 


2.1 


X 1.6 - 


-14° 


Mkn 231 


25 


1.0 


228 


1.8 


X 1.6 


9° 


Mkn 348 


24 


1.3 


287 


1.8 


X 1.5 


47° 


Mkn 668 


15 


2.0 


736 


1.9 


X 1.7 - 


-53° 


NGC 1068 


100 


0.9 


1395 


2.8 


X 1.6 - 


-23° 


NGC 2110 


114 


0.9 


126 


2.7 


X 1.5 - 


-35° 


NGC 2639 


60 


0.8 


85 


1.8 


X 1.7 


68° 


NGC 2992 


90 


0.7 


49 


3.7 


X 1.5 


17° 


NGC 3079 


141 


0.6 


83 


1.7 


X 1.6 


21° 


NGC 3504 


29 


1.4 


52 


1.8 


X 1.7 


75° 


NGC 4151 


49 


0.8 


193 


1.9 


X 1.8 


70° 


NGC 5506 


89 


0.7 


223 


2.0 


X 1.5 


9° 



Note. — The reported continuum measurement is the peak flux density 
of the continuum image but not necessarily the integrated continuum 
emission from the source. 



-44- 



Table 3. Hi Absorption Properties 



Source Peak r J S^, dv / dv Nhi Separation Proj. Sep. 

(Jy km s'l) (km s'^) (lO^i cm'^) (") (pc) 

(1) (2) (3) (4) (5) (6) (7) 



Mkn 3 


< 0.005 




< 0.09 




< 0.10 




< 0.02 










Mkn 6*^ 


0.45± 0.01 


1.1± 


0.06 


14.0± 


1.6 


2.6± 


0.3 


l.OiO.lO 


370 ± 


36 


Mkn 231" 


0.073 


0.004 


3.6 


0.17 


33.9 


3.8 


6.3 


0.7 


0.10^0.44 


82-^ 


360 


Mkn 348 


< 0.004 




< 0.08 




< 0.27 




< 0.05 










Mkn 668 


< 0.006 




< 0.12 




< 0.16 




< 0.03 










NGC 1068 


0.078 


0.004 


1.5^ 


2.9 


2.9^ 


15 


0.53^ 


2.73 


1.6^ 8.4 


120^ 


625 


NGC 2110 


0.120 


0.012 


1.0± 


0.16 


17.3± 


2.1 


2.5± 


0.4 


1.5±0.16 


230 ± 


24 


NGC 2639 


< 0.003 




< 0.05 




< 0.56 




< 0.10 










NGC 2992 


0.212 


0.020 


1.1 


0.12 


31.2 


3.4 


5.7 


0.6 


1.0 0.13 


152 


19 


NGC 3079 


0.788 


0.007 


9.3 


0.19 


135.3 


2.3 


24.6 


0.5 


<0.1 •■• 


< 3 




NGC 3504 


0.188 


0.026 


1.3 


0.16 


30.2 


3.0 


5.5 


0.7 


<0.1 •■■ 


< 11 




NGC 4151^ 


0.241 


0.024 


0.9 


0.05 


21.0 


3.0 


3.8 


0.6 


<0.1 •■• 


< 6 




NGC 5506 


0.109 


0.003 


3.0 


0.10 


15.1 


0.5 


2.6 


0.1 


<0.1 •■■ 


< 11 





Note. — Columns: (1) Source ID; (2) peak 21 cm opacity; (3) absorption line strength; (4) the 
velocity-integrated opacity of the strongest absorption line; (5) the corresponding H I column density, 
assuming T5 = 100 K; (6) the angular separation between the AGN and the strongest Hi absorption; 
(7) the projected separation in parsecs. Note that the measurements of r are necessarily lower limits 
because of confusion and unknown covering fraction. Ranges (^) are given for the cases where the Hi 
absorption is spatially resolved. For sources with unresolved H I absorption or continuum, the equivalent 
width and column density are lower limits. Expressed limits are 3o", assuming a linewidth matching the 
~ 20 km s~^ channel-width. For the NGC 4151 upper limits on Hi separation, we assumed that the 
AGN is located at C4W, the central VLBI component, after Mundell et al. (1995). 

'^Values include MERLIN measurements from [Mundell et al. (1995|) . 
^Values include MERLIN measurements from [Gallimore et al. (19981 ). 
'^Values include VLBA measurements from [Carilli, Wrobel, fc Ulvestad (19971). 
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Table 4. H l Absorption Line Profiles 



Source 


^max 




FWHM 


^^50 




^50 - V. 


max 


'^rnax 


^sys 


^50 - Vsys 




(km S~ 




(km s 




(km s~ 




(km s" 




(km s 




(km 




(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




Mkn 6 


5572.2±0.6 


57.7±4.2 


5575.1±4.2 


-2.9±2.6 


-67.8±10.0 


-64.9±4.9 


Mkn 231 


12664.2 


5.9 


211.8 


6.8 


12647.5 


6.8 


16.7 


6.5 


14.2 


20.9 


-2.5 


9.4 


NGC 1068 


1117.8 


1.5 


161.1 


5.0 


1170.2 


5.0 


52.4 


4.7 


-30.2 


5.2 


22.2 


6.8 


NGC 2110 


2626.3 


4.3 


88.5 


7.6 


2616.3 


7.6 


10.0 


5.6 


291.3 


20.5 


281.3 


9.4 


NGC 2992 


2296.5 


0.9 


95.0 


4.8 


2326.6 


4.8 


-30.1 


4.5 


-17.5 


6.1 


12.6 


6.5 


NGC 3079 


1136.2 


0.4 


132.5 


1.4 


1130.0 


1.4 


6.2 


1.2 


12.2 


10.0 


6.0 


1.8 


NGC 3504 


1542.0 


1.0 


135.5 


7.4 


1525.3 


7.4 


16.7 


6.0 


8.0 


2.2 


-8.7 


9.5 


NGC 4151 


994.7 


0.9 


89.8 


5.2 


991.0 


5.2 


3.7 


3.3 


-2.3 


3.1 


-6.0 


6.1 


NGC 5506 


1802.3 


0.8 


95.0 


3.5 


1821.7 


3.5 


-19.4 


3.5 


-12.7 


9.0 


6.7 


5.0 



Note. — Columns: (1) Source ID; (2) velocity at absorption line minimum (i.e., peak line 
strength); (3) full-width at half maximum; (4) the velocity mean across the FWHM; (5) as a 
line asymmetry parameter, the difference between v^ax and ^59, with positive values indicating 
a red asymmetry in the line wings; (6) & (7) the difference between the systemic velocity and 
the line profile velocities, with positive values indicating a red-shift relative to systemic. 
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Table 5. Soft X-ray Absorbing Columns 



Source 


Nh 

(1021 em-2) 


Reference 


Mrk 3 


707.9tt^ 


(1) 


Mrk 6 


38.01^:° 


(2) 


iV±l K Z O X 


c q+3.3 




Mrk 348 


112.21^0 


(1) 


Mrk 501 


< 18. 


(4) 


NGC 1068 


> 1000. 


(5) 


NGC 2110 


28.4+}-^ 


(6) 


NGC 2992 


1-31^:? 


(6) 


NGC 3079 




(7) 


NGC 3227 




(8) 


NGC 4151 


22.61}° 


(8) 


NGC 5506 


21.911^ 


(9) 



Note. — For those cases where 
more than one X-ray spectral model is 
available in the literature, we selected 
the model giving a reasonable fit to 
the data and the best constraint on 
the hydrogen column. References: 
(1) Mulchaey, Mushotzky, & Weaver 
(1992); (2) Feldmeier et al. (1999); 
(3) Wang et al. (1996); (4) Mushotzky 
et al. (1978); (5) Matt et al. (1997) 
and references therein; (6) Turner et al. 
(1997); (7) Read, Ponman, & Strickland 
(1997); (8) George et al. (1998); (9) 
Mulchaey et al. (1993). 
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Table 6. Model Disk Parameters 



Source 


i 


PA 


Vrot 


A?"/// (face-on) 




^out 








(deg) 


(deg) 


(km s~^) 


(xlO^i 


cm ^) 


(pc) 


(pc) 


(M©) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Mkn 6 


51 


130 


140 ±4 


1.65 


±0.07 


293 •• • 


426 •■ • 


6.60 


143.8 


NGC 1068 


40 


-45 


320 82 


0.82 


0.056 


212 •• • 


957 •■ • 


7.26 


64.5 


NGC 2110 


53 


161 


412 5 


2.31 


0.25 


250 ■• ■ 


552 •• • 


7.15 


57.0 


NGC 2992 


70 


-154 


147 8 


0.73 


0.05 


36 ±34 


921 ■• ■ 


7.20 


185.9 


NGC 3079 


84 


166 


177 1 


1.53 


0.01 


8 1 


282 ±4 


6.48 


1089.0 


NGC 3504 


22 


-33 


235 27 


2.33 


0.21 


8 1 


< 312 ■•■ 


< 6.76 


50.5 


NGC 4151^ 


21 


22 


-148 8 


9.64 


0.63 




7 ... 


4.85 


157.8 


NGC 4151'^ 


50 


180 


839122 


10.20 


0.69 




10 •■• 


4.87 


66.6 


NGC 5506 


80 


-87 


181 5 


0.25 


0.01 


7 2 


406 26 


6.00 


963.2 



Note. — Columns are: (1) Source ID; (2) inclination; (3) position angle to receding major 
axis; (4) the deprojected rotational velocity at 100 pc radius; (5) face— on column density through 
the disk; (6) the deprojected inner and (7) outer radii; (8) logarithm of the Hi mass (in Mq) 
assuming a uniform and symmetric disk; and (9) x^- Error-bars are given for fitted parameters. 
Where not explicitly fit, the inner and outer radii span the deprojected range of the observed 
absorption. The reported values are for 50 spectral data points and 3-5 fitted parameters. We 
fit the MERLIN data of Gallimore et al. (1998) for Mkn 6; the values are for 53 spectral data 
points and three fitted parameters. Note that the Hi masses of the model disks are less than 
0.2% of the dynamical mass at the corresponding radii. References for the kinematic axes are, 
Mkn 6: Meaburn, Whitehead, & Pedlar (1989); NGC 1068: Brinks et al. (1997) and Galhmore 
et al. (1994); NGC 2110: Wilson & Baldwin (1985); NGC 2992: Colina et al. (1987); NGC 3079: 
Irwin & Seaquist (1991); NGC 3504: Kenney, Carlstrom, & Young (1993); NGC 4151: Pedlar et 
al. (1992) (galaxy disk) and Mundell et al. (1995) (AGN disk); NGC 5506: Maiohno et al. (1994). 

^The inclination and position angle are fixed to those of the outer galaxy disk. 

^The inclination and position angle are fixed to those predicted for an AGN disk (i.e., 
perpendicular to the radio jet). 
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Table 7. Streaming Velocity Models 



Source 






M 




Ax' 




(km s" 




(M^yr-i) 






(1) 


(2) 




(3) 


(4) 


(5) 


Mkn 6 


< 59.0 




< 0.57 ■ • ■ 


95.4 


48.4 


NGC 1068 


-46.6±11.2 


-0.90±0.21 


36.8 


27.7 


NGC 2110 


47.3 


11.7 


1.16 0.22 


42.4 


14.6 


NGC 2992 


37.6 


5.4 


0.67 0.09 


131.8 


54.1 


NGC 3079 


16.6 


0.6 


0.18 0.01 


927.1 


161.9 


NGC 3504 


9.6 


19.4 


0.18 0.36 


50.5 


0.3 


NGC 415P 


-0.32 


0.04 


-182.4 22.5 


48.5 


109.3 


NGC 4151^ 


-0.07 


0.01 


-24.8 3.8 


46.1 


20.5 


NGC 5506 


19.3 


0.9 


0.05 0.01 


715.8 


247.4 



Note. — Streaming velocities based on disk model 
fitting to the H I absorption profiles and locations. Positive 
velocities imply motion toward the center of the galaxy. All 
models also include a rotational component, but the fitted 
values are not significantly different from those listed in 
Table ^. The exception is Mkn 6, for which the fitted 
rotational velocity drops to ~ 40 km s~^; the fitted value 
therefore places a limit appropriate for gas dominated 
by radial motion. The columns are (1) Source ID; (2) 
streaming velocity; (3) corresponding mass infall rate; (4) 
x'; and (5) Ax', the improvement in x' relative to the 
purely rotational model (Table ||). 

^The inclination and position angle are fixed to those of 
the outer galaxy disk. 

''The inclination and position angle are fixed to those 
predicted for an AGN disk (i.e., perpendicular to the radio 
jet). 
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